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Broad redshifted line as a signature of outflow 

Lev Titarchuk^'^, Demos Kazanas^ and Peter A. Becker ^ 

ABSTRACT 

We formulate and solve the diffusion problem of line photon propagation 
in a bulk outflow from a compact object (black hole or neutron star) using a 
generic assumption regarding the distribution of line photons within the outflow. 
Thomson scattering of the line photons within the expanding flow leads to a 
decrease of their energy which is of first order in v/c, where v is the outflow 
velocity and c is the speed of light. We demonstrate that the emergent line 
profile is closely related to the time distribution of photons diffusing through the 
flow (the light curve) and consists of a broad redshifted feature. We analyzed the 
line profiles for the general case of outflow density distribution. We emphasize 
that the redshifted lines are intrinsic properties of the powerful outflow that are 
supposed to be in many compact objects. 

Subject headings: accretion — stars:radiation mechanisms:nonthermal 

1. Introduction 

The problem of photon propagation in an optically thick fluid in bulk motion has been 
studied in detail in a number of papers. The idea that photons may change their energy 
in repeated scatterings with cold electrons in a moving fluid was suggested more than 20 
years ago by Payne & Blandford (1981), hereafter PB81 and Cowsik & Lee (1982). This 
process, often referred to as dynamical or bulk flow Comptonization, is similar in many ways 
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to Comptonization by hot electrons once the thermal velocity is replaced by bulk velocity v; 
there is however a qualitative difference, in that their energy gain is linear in the velocity v, 
rather than quadratic as is in the case of Compton scattering in a medium at rest. Photons 
diffusing and scattering in a medium with bulk flow can gain or lose energy to the flow 
depending on the divergence of its velocity field: they gain energy for V ■ v < 0, while they 
lose energy in the opposite case. The case of a converging radial infiow, where V • v < 0, 
has been treated by (among others) PB81, Nobili, TuroUa, & Zampieri 1993, TuroUa et 
al. (1996) and Titarchuk, Mastichiadis & Kylafis (1996; 1997), (hereafter TMK96, TMK97 
respectively). These investigations showed that, when monochromatic radiation with v — 
is injected at large Thomson depth in a spherical inflow, the emergent spectrum develops 
a broad, power-law tail at z/ > z/*. The power law index is related to a combination of the 
flow Thomson depth and its velocity gradient, and this has been invoked to explain the high 
energy spectra of accreting BH candidate sources in their "high" state (Laurent & Titarchuk 
1999; Shrader & Titarchuk 1999). 

In this Paper we treat the problem of an outflow and show that a similar broad spectrum 
is formed but now at energies v < v^,, with the power law index dependent on the velocity 
gradient. This problem is also relevant observationally because of the well known presence 
of red wings in the fluorescent Fe Ka line proflles observed in the AGN and galactic BH 
spectra. TuroUa et al. (1996) were first who analyzed the radiative transfer in an expanding 
atmosphere. They found a particular solution of the PB81 equation for diverging fiow using 
Fourier transformation. They showed that the overall behavior of the spectra is similar 
to that of the converging fiow but somehow reversed, since now photons can drift only to 
frequencies lower than v^,. Furthermore, they found that adiabatic expansion produces a 
drift of injected monochromatic photons towards lower frequencies and the formation of a 
power-law, low energy tail with spectral index 3. 

Photons can also lose energy in scattering with cold electrons due to recoil. This process 
was treated by Sunyaev & Titarchuk (1980; hereafter ST80), who demonstrated that the 
emergent spectrum of the high energy photons injected in the static cold electron cloud is 
determined by the initial photon distribution in radius, which in turn dictates the photon 
escape probability distribution in time. The resulting spectrum for monochromatic injection 
is a redshifted feature of spectral width Ai/ which depends on the line frequency and the 
Thomson depth tq, i.e. Ai/ ~ — i/*(/iz/*/meC^)rQ . The observation of red-skewed Ka— lines 
observed in a number of AGN (and galactic black hole systems) (Fabian et al. 1995), led to 
the consideration of the role of this process in producing the observed Fe fine profiles. This 
was however dismissed by Fabian et al. (1995), who correctly pointed out that the required 
Thomson downscattering of the line photons would also introduce a (non-observed) break in 
the spectrum at energy E ~ rrieC^ /tq ~ 20 keV. An approach along the same lines by Misra 
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Sz Kembhavi (1998) was rebutted on the basis of line and continuum variability observations 
by Reynolds & Wilms (2000). 

Following these works, the prevailing process for producing the observed profile of the 
Fe Ka— line was considered to be Doppler broadening by the kinematics of a cold {kT < 0.25 
keV) accretion disk reprocessing the X-ray radiation of an overlying corona into fluorescent 
Fe KcK— line and reflection continuum. It was shown (Fabian et al. 1989) that such an 
arrangment could produce lines of the desired width and asymmetry provided that the 
disk extended to its innermost stable orbit (3 Schwarzschild radii). However plausible, this 
interpretation is not entirely without problems; for instance, the (occasionaly) exceedingly 
broad red wing of the line observed requires (in such cases) that the disk illumination be 
concentrated very close to its inner edge {F^ oc r"^; Nandra et al. 1999) more than most 
models would allow; this led Reynolds & Begelman (1997) to consider fluorescent emission 
from within the innermost stable orbit, by matter inspiraling into the black hole. Also, as 
indicated by the models of Nayakshin et al. (2000); Ballantyne et al. (2001), the ionization 
of such a disk by the intense X-ray radiation might invalidate some of basic assumptions 
associated with this interpretation. Of additional interest is also the fact that there is a 
marked absence of a blue-shifted wing in the line, a feature expected for a random orientation 
of accretion disks and the observers' lines of sight. 

Motivated by the above facts and the importance of the Fe Ka— line in probing the strong 
field limit gravity we believe it is important that all potential alternatives to line broadening 
be considered in detail. To this end we examine in the present Paper under what conditions 
these broad Fe line profiles could be also attributed to the effects of an outfiow rather than 
solely to accretion disk kinematics. The formulation of the photon diffusion problem in an 
outflow and its solutions are described in §2, and in §3 our results are discussed and certain 
conclusions are drawn. 



Let N{r) — No{ro/r)^ be the radial number density proflle of an outflow and let its 
radial outward speed be 



2. 



Radiative Transfer in a Bulk Outfiow 



2.1. Main equations and parameters 



Vb/c={Mout/'^7rcNorl){ro/r) 



2-/3 _ 



mout{ro/r) 



2-/3 



(1) 
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obtained from mass conservation in a spherical geometry (here Mout — 'inr'^VbN) . The 
Thomson optical depth of the flow from some radius r to infinity is given by 

/oo 
N,{r)aT dr = aTNoro{ro/rf-' / {P - 1), (2) 

where Ne{r) = N{r) is the electron density, is the Thomson cross section, tq is a radius 
at the base of the outflow (this deflnition of optical depth makes the tacit assumption that 
P>1) and tt,o - T(ro) - arNoro/iP - 1). 

The transfer of radiation within the flow in space and energy is governed by the photon 
kinetic equation (BP81, Eq. 18) for the photon occupation number n(r, i/), which in steady 
state reads 

1 c 1 Ofi ~ 

-Vb • Vn + -V • (-Vn) + 3 (V • Vb) z/— = -j{r, v), (3) 

where k = Ne{r)(TT is the inverse of the scattering mean free path, Vb = Vi,er, is the flow 
velocity, Br is the radial unit vector and j(r, u) is the photon source term. 

The spectral flux F(r, u) (PB81) is given in terms of n(r, u) by 

1 1 Oti 

and must satisfy the following boundary conditions: (a) Conservation of the frequency inte- 
grated flux over the outer boundary, namely 

J F{r, v)dv oc as r — > 00. (5) 

(b) The occupation number should be zero at the inner boundary r — Tq, where the bulk 
outflow starts and the outflow radiation is zero, i.e. 

n(ro) = 0. (6) 



2.2. Solution for a general case of any optical depth tq 

According to a theorem (Titarchuk 1994, appendix A), the solution of any equation 
whose LHS operator acting on the unknown function n(r, v) is the sum of two operators 
and Ly, which depend correspondingly only on space and energy and the RHS, j{r,i'), is 
factorizable, i.e. 

LrU + L^n = -j(r, ly) = -f{r)if{p) (7) 
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with boundary conditions independent of the energy i.e. 

L^^^n = as r — >• oo, Lf^'n = for r = ro , (8) 

is given by the convolution of the solutions of the time-dependent problem of each operator, 
namely 

poo 

n{r,iy)^ P{r,u)X{iy,u)du. (9) 
Jo 

Above u = NoCT^ct is the dimensionless time and P(r, u) is the solution of the initial value 
problem of the spatial operator Lr 

^-L,P, P(r,0) = /(r) (10) 

with boundary conditions 

L^p^O as r^oo, L^^^P = at r = ro, (11) 

and Xlu, u) the solution of the initial value problem of the energy operator 

dX 

— = L.X, X(z/,0) = ^(z/) (12) 

with boundary conditions 

v^X when i/ ^ 0, oo. (13) 

In the case where the photon energy change is due to the electron recoil (see details in 
ST80), we have 

-Vn 



L^n = — V 
3k 



K 



(14) 



and 



1 a . 4 



where z = hu/meC^. The solution of the time- dependent problem for the energy operator 
Li, (the Green's function) and ip{z) = S{z — z^)/ z^ is (see ST80) 

u \ ivi \ { 6[{l/z-u)~^ - z^]/z if 2; < 2;* and < M < 1/2; , . 

Jiz.u) — zXiz.u) — { ' , . 16 

^ ' ^ ^ ' ^ \ otherwise, ^ ^ 

where J{z, u) — z^X{z, u) is the radiation intensity. Substitution of Eq.(16) into Eq.(9) gives 
the form of the response function (the Green's function) to monochromatic line injection in 
a static cloud _^ 

G,S^) = z^n{rouu z) = —P[rouu {^/z - l/z.)]. (17) 
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For a local photon injection the time-dependent diffusion distribution in a bounded medium 
has the form of a fast rise and exponential decay (see ST80 and Wood et al. 2001), where 
the characteristic dimensionless time {u = aTNoct) for the exponential decay is always 
proportional to t^q, independent of the source distribution. Thus the asymptotic form 
of the redshifted wing of the downscattering line (i.e. for u <^ z/*) is always exponential, 
Gz,{z) oc exp{—Ao/zT^ q). It is evident from this formula that a significant line redshift 
can occur only when hp/rrieC^T^ o ~ where Aq is of order of a few. 

To obtain the Green's function of the same problem for a diverging flow we use the 
same method applied in the downscattering case. The time-dependent problem for the 
energy Green's function is (see Eq.3 and Eqs. 12-13) 

X{iy, 0) = S{iy - v,)lv^ (19) 

with boundary conditions 

v^X when i/ ^ 0, oo. (20) 

The problem (18-20) is an initial value problem for the first order partial differential equation 
and it can be found using the method of characteristics (see e.g. Zel'dovich & Raizer, 2002) 
Because the occupation number X{y^ 0) is conserved along the characteristics u — hi.VQlvoi 
equation (18) the solution of the problem (18-20) is 

J^X^.n) = i/2X(i/,M) = e-2"5(t/e" - t/,). (21) 

Substitution of J^^{u,u) from Eq.(21) into Eq.(9) gives us the Green's function 



G,X^) = i^Mrout, 1^) = - (-) P[rout, H^*/i^)]. (22) 



The time distrubition P{rout, tj) can be found from the solution of the time-dependent prob- 
lem of the space operator and initial condition (see Eq. 10) 

du ' V • (Vb/c) V • (Vb/c) ' ^ ^ 

P(r,0) = /(r) (24) 

with two boundary conditions 

-{1/k)VP (xr~^ at r^oo, P = at r = ro. (25) 
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Introducing r = tt as a spatial variable Eq.(23) can be rewritten as follows 

"a^P 2 dP 



dr^ {(5 - l)r dr 



dP 

+ 3r-^ (26) 



where = [Mout/^T^cN^rDril^'^^"-^-'^ = moutri!'-^^'^^~'\ 1 = {W - 3)/(/5 - 1). In order 
to solve the above time-dependent problem we look for solution of the eigenvalue problem 
for the operator Lj. (Eq. 26) with the appropriate boundary conditions, namely 

LrR + X^R = 0, (27) 

with the boundary conditions 

dR/dr oc r^/^^-^) oc r'^ as r ^ (or r ^ oo) and P = at r = Tq. (28) 

The nontrivial solution of equation (27) satisfying the boundary conditions (28) (i.e. eigen- 
fuction) is 

i?,(r) - r^P+^)liP-^)m3 + 1 + /3 + 2, -3a,(/3 - 1)tV(/5-i)], (29) 

where $(a, 6, 2;) is the confluent (or degenerate) hypergeometric function. The boundary 
conditions (28), along with formula (29), implies that the eigenvalues are roots of the equation 

$(/3 + 1 + A2/3/3, /3 + 2, -Te//,o) = 0, (30) 

where Te//,o = 3moMt(/3 — 1)to- Using the asymptotic form of the confluent hypergeometric 
function $(0,6,2;) for large arguments 1 (Abramowitz & Stegun 1970, Eq.[13.1.4]) 

and equation (30) we can find the eigenvalues for Te//,o = 3(/3 — l)mo„tro ^ 1 as roots of the 
following equation 

l/r(l- A^/3/3) =0, (32) 

where r{z) is the gamma-function. Because r{z) goes to 00 for nonpositive integers, z — 
— {k — 1), {k — 1, 2, 3,...), the eigenvalues are 

= 3k/ p. (33) 

If the monochromatic line sources are distributed according to the eigenfunction /(r) = 
Rkir) then P{t,u) is 

Pfe(T,ii) = Pfe(T)exp(-A^ii). (34) 
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Thus when /(r) oc -Rjt(r) the Green's function Gv^{v) [see formula (22)] is a power law 



1 



G.X^)^-\-\ . (35) 
For Te//,o ^ 1 the spectral index 

= 3 + = 3(1 + A;//3) (36) 

and the least spectral index q; = q;i = 3(1 + 1//5). In Figure 1 we present numerical 
calculations of the least photon index T — a — 1 using Eq. (30) for various values of (3 and 
Te//,o- As it is seen from this figure F— values converge to the asymptotic values Fas = 2 + 3//3 
for Te//,o > 1. 

In the general case of a given source distribution /(r), rriouti ^o, and /? one should 
calculate the eigenvalues using Eq. (30) and the expansion coefficients bk using the following 
formulae 

r^lj{r)Rkir)fir)dr/Hkiro), (37) 

Hk^ r ^P{T)Rl{T)dT, (38) 
Jo 

where V'(t) = r^^~'^l^y exp(^) is a weight function, ^ = 3C^^(/3 — I'^r^/^/^'^^ is inversely 
proportional to the radius r (see Eq. 2) and Cy^ = mout^o^'"^^^^^'^^ ■ It is worth noting that 
the variable ^ is equal to the trapping radius for advection, rt, divided by the radius r, where 
rt — {GM / c^){Moutl Me), with Me — A^T^GMmn / cut denoting the Eddington accretion rate 
(see e.g. Becker & Begelman 1986, TMK97). 

Equations (37-38) are obtained using the orthogonality properties of the eigenfunctions 
i?fe(r). In appendix A we show that the square norm can be computed analytically in 
closed form, namely that 

where pij) = [(/? — l)//3Ct,Jr'''-(/'(r) [see definition of 7 in Eq. (26) above]. Convenient 
formulae for calculations of derivatives of Rk above are given in Appendix A. The numerator 
of equation (37) Ck [a scalar product of /(r) and Rk{T)] can be also calculated analytically 
in the case of the 5— function source distribution [/(r) = S{t — t^), where < Tq]: = 

ip{n)Rk{n). 

The formula for the flux J-'{i') = [(r/ro)^F(r, i/)] |r^oo (see Eq. 4) is 
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In order to calculate Rk for A; » 1 in equations (29, 37-38) one should use the asymptotic 
form of the confluent functions $(a, b, z) for large values of a and z > 1 (Abramowitz & 
Stegun 1970, Eq. [13.5. 13]). 

$(a, 6, z) = r(6)e"(Z/4)(V4-V2)7r-V2 ^os (Z^^ _ 67r/2 + 7r/4), (41) 

where Z = (4a — 2b) {—z). In Figure 2 we present the results of calculation of the time 
distribution V{u) oc [r^/(^~'^^(?P/c?r[r, |^^o using the expression [see Eq. (34) for the 
definition of Pk{T, u)] 

oo 

P{T,u) = J2hPk{r,u) (42) 

k=l 

and equations (29-30, 37-38). These particular calculations are made for a constant velocity 
wind (/? = 2), TJiout — 0.8 and tq = 2, using as initial conditions a monochromatic source 
[9?(z/) = (5(z/ — z/*)] and four different spatial distributions given by /(r) = exp {rj^), with 
r] — 0, 1.5, 2.5 (i.e. sources increasingly concentrated towards the base of the flow for 
increasing 77) and /(r) = Ri{t), i.e. a distribution according to the first eigenf unction. 

All spectra are related to P{t,u) through the transformation of Eq. (22). In Figure 3 
we show the evolution of the emergent spectra as a function of the initial source distribution 
/(r). These calculations are made using expression (40). 



2.3. Solution for the asymptotic case tq » 1 

TuroUa et al. (1996) hereafter TZZN, analyzed this case in detail and they found 
analytical solution of this problem using Fourier transformation of the main equation (3). 
Below we show the asymptotic of our solution for tq ^ 1. 

Because f3\k/3 = k for r ^ 1 (see formula 33) we can rewrite formula (29) for the 
eigenfunction as follows 

i?,(r) = + i + k,P + 2,-^]. (43) 

The confluent hypergeometric function $(a, b, x) assumes the transformation [see Grad- 
shteyn & Ryzhik 2000 (hereafter GROO), equation 9.212.1]: 

$(6-a,6, -x) = e-^$(a,6,x). (44) 

In our case a = 1 — k and b = (3 + 2. Then we can express Rk through Laguerre polynomial 

i?.(r)=r(^+^)/(^-^)e-^LK(0 (45) 
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keeping in mind that 



$[l-fc,/3 + 2,e]cxLK(0 



(46) 



(see GROO, equation 8.972.1). We neglect the constant factor in front of the right hand side 
of relation (46) to derive formula (45) because the eigenfunctions Rk are always determined 
within a constant factor. 

Thus the square norm of Rk is: 

/■oo fOO 

Hu^ / e-V/(^-^)[LK(0]^dr = (/3-l)S-(/^+^) / e-^e^+^[LK(0]^de 
Jo Jo 

We transform the integral above using the definition of ^ = ^Cvi,{P — 1)t^/^'^~^^ = ^t-i/C/^-i) 
and we obtain 

//,^(/3-1)E-(/^+^) ^^^+_\|,^\ (47) 

where we have used equation (7.414.3) from GROO, to evaluate the integral above for case 
n — m^. It is evident that the constant factor in the right hand side of equation (47) is 
independent of k and depends on (5 and C^^ only. 

In order to determine the expansion coefficient 6„ for the 5— function source distribution 
[/(t) = 5{t — T*)] we need to integrate the numerator in Eq. (37) in r over a small region 
around t — < Tq. Taking into account Eq. (47) this yields 



6feOc(^-l)!/r(/3+l + A;)Lf+}(e.) 



(48) 



where = 3C^,(/3 - 1)tI'^'^~^^ and we have neglected the constant factor in the right hand 
side of equation (48) which is independent of k and depends on (5 and r* only. Our solution 
for the flux T is therefore given by (see Eq. 40) 



1 / \ 3 oo 



{k-l)\ 



(49) 



In fact, using formula (8.976.1) of GROO, for this series the final, closed-formed solution for 
the Green's function can be written as 



jr(i/) oc — — 



3(l+l//3) 



exp 



1 - (i^M)^/'^, 



1-1^ 



3//3' 



-(/3+2) 



(50) 



*It is worth noting that in the fifth edition of Gradshteyn & Ryzhik' table of the integrals there is a typo 
in formula (7.414.3). 
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where in formula (8.976.1) of GROO we have taken the hmit y — > and used the asymptotic 
form of the modified Bessel function = {w /2Y^^ /V[j3 + 2) for the argument w <^1. 

The shape the spectrum (Eq. 50) is not identical to that obtained by TZZN. In fact, the 
low-frequency limit of the spectrum is a power law 



with the index a — ai — 3(l-|-l//3). The difference between our result and that of TZZN can 
be explained by the different boundary conditions in two cases. It is easy to show that the 
TZZN solution obtained using the Fourier transform technique is a particular solution of Eq. 
(3) with a reflection inner boundary condition, while our solution is related to an absorption 
boundary condition. A reflection boundary condition, by conserving the number of photons, 
as docs the pure scattering involved in the process we examine, leads, not surprizingly, to a 
Wien-like spectrum, i.e. to a = 3. The situation is different with an absorptive boundary 
which does not conserve the number of photons, leading to a harder spectrum as the boundary 
is more likely to absorb the photons that scatter longer, i.e. the lower energy ones. It is 
therefore not by chance that the TZZN index is 3 but our index is 3(1 + 1//3). A similar 
situation is obtaind for the spectra of converging flows. The spectral indices are different 
depending on cither reflection or absorption inner boundary conditions assumed. TMK96 
and TMK97 (see also Turolla, Zane & Titarchuk 2002) flnd that the spectral indices of the 
power-law are close to and 2 in the reflection and absorption cases respectively. 



We have presented above an idealized treatment of the radiative transfer of monochro- 
matic photons within a divergent flow, a problem of also observational interest, in view of 
the observations of broad, redshifted Fe lines in the extragalactic and galactic black hole 
candidate spectra. Not surprisingly, the diverging flow leads to a redshift of the line the 
photon energy to produce lines with broad, red wings, not unlike those of the Fe Ka lines 
observed. In our view, one of the more interesting aspects of our proposal and the results of 
our calculations is the natural absence of a corresponding blue wing in this feature without 
the need to invoke a specific geometrical arrangement for the emission. 

Jets and outflows as a means for producing the observed broad red wings of the Fe Ka 
fines have also been considered by Fabian et al. (1995). However, these authors had a rather 
different mechanism in mind: excitation of the line in an optically thin, cold jet by external 
X-ray illumination. The redshifted wing in this case would, they correctly argued, ought to 
be effected by the transverse Doppler effect or a gravitational red shift. We concur with their 




(51) 



3. 



Discussion and Conclusions 
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assesment that such and set-up for producing the observed hne profiles would require a rather 
special geometric arrangement, as even a large (but smaller than 90°) angle of the outflow to 
the observer's line of sight should result in blue shifted lines which are generally not seen (see 
however Pounds et al. 2001). Our proposed process is different in that the line is produced 
in an (effectively) optically thick medium. Its red wing is the result of multiple scattering 
and the associated first order redshift in each such scattering. This process produces a red 
wing to the line without a particularly fine tuned geometric arrangement. 

As indicated by our calculations, the most crucial parameter in determining the line 
shape is the effective scattering depth of the flow Te//,o and the photon source distribution. 
Clearly, a very extended source distribution allows for a large number of photons to be 
produced at regions of very low scattering depth (case of ?] = in Figure 3) leading to a 
rather narrow line profile. An increase in the value of rj (i.e. producing a larger fraction of the 
photons deep into the flow) leads to much broader profiles which can even exhibit a shoulder 
or a secondary minimum, not unlike the Fe Ka line feature that was interpreted as evidence 
for infalling matter in NGC 3516 (Nandra et al. 1999). The origin of this feature is well 
understood within our model: the line narrow core corresponds to photons produced at low 
values of ^ (large distances r), which escape without much scattering in the flow, while the 
broader, redshifted profile is formed by the photons (trapped in the fiow) produced at high 
values of ^ which diffuse their way out through the expanding flow. This behavior can be 
guessed from Fig. 2 where the time distribution of photons escaping the flow is presented: 
There is a narrow peak near zero, representing the photons escaping directly (essentially 
with little energy change) and a broader peak comprising the photons diffusing through the 
flow while at the same time suffering the energy losses that distribute them on the red wing 
of the line. It is worthwhile emphasizing that the narrow core would completely vanish in 
the spectrum for the asymptotic case of tq ^ 1 and only the broad feature followed by the 
low frequency power-law arc present there. Wc should also point out that photon scattering 
in diverging flows introduces besides redshifted line proflles also negative time lags in the 
spectra (see Figs. 2-3 and formula 22). In fact, a photon of initial energy i/* loses energy on 
the way out. If it escapes after time u {u scatterings) its flnal energy will he u — e"'^!/^ and 
it will escape later (how much that depends on the speciflcs of the scattering plasma) than 
a photon of energy i/* which has not scattered at all. We suggest the combination of these 
features to be intrinsic signatures of any diverging flow. 

Finally, while our results are for simplicity speciflc to a particular flow proflle, different 
density and velocity proflles associated with more general types of outflow, such as ADIOS 
(Blandford & Belgelman 1999), could in fact offer much broader parameter space to explore 
with more diverse results. We defer such an exploration to a future publication. 



-13- 



L.T. acknowledges the support of this work by the Center for Earth Observing an 

Space Research of George Mason University. D.K. would like to acknowledge support by a 
Chandra GO grant and would like to thank Frank C. Jones for helpful discussions. We also 
acknowledge the thorough analysis of this paper by the referee and his/her constructive and 
interesting suggestions. 

REFERENCES 

Abramowitz, M. & Stegun, L. 1970, Handbook of Mathematical Functions (New York: 
Dover) 

Ballantyne, D.R., Ross, R.R. & Fabian, A.C. 2001, MNRAS, 327, 10 
Becker, R A., & Begelman, M. C., 1986, ApJ, 310, 534. 
Blandford, R.D. & Payne, D.G. 1981, MNRAS, 194, 1033 
Blandford, R.D. & Begelman, M.C. 1999, MNRAS, 303, LI 
Cowsik, R. & Lee, M.A. 1982, Proc.R.Soc. London A, 383, 409 
Fabian, A.C. et al. 1995, MNRAS, 242, 14P 
Fabian, A.C. et al. 1989, MNRAS, 277, Lll 

Gradshteyn, I.S. & Ryzhik, I.M. 2000, Table of Integrals, Series and Products (sixth edition; 
Eds. A. Jeffrey and Z. ZwiUinger) (New York: Academic Press) (GROO) 

Laurent, P. & Titarchuk, L. 1999, ApJ, 511, 289 

Misra, R. & Kembhavi, A.K. 1998, ApJ, 499, 205 

Nandra, K. et al. 1999, ApJ, 523, L17 

Nayakshin, S., Kazanas, D. & Kallman, T.R. 2000, ApJ, 537, 833 
NobiUi, L., TuroUa, R. & Zampieri, L. 1993, ApJ, 404, 686 
Payne, D.G. & Blandford, R.D. 1981, MNRAS, 196, 781 (PB81) 
Pounds, K. et al. 2001, ApJ, 559, 181 
Reynolds, C.S. & Begelman, M.C. 1997, ApJ, 488, 109 



-14- 



Reynolds, C.S. & Wilms, J. 2000, ApJ, 533, 831 
Shrader, C.R. & Titarchuk, L. 1999, ApJ, 521, L121 
Sunyaev, R.A. & Titarchuk, L.G. 1980, A&A, 86, 121 (ST80) 
Titarchuk, L. 1994, ApJ, 434, 570 

Titarchuk, L., Mastichiadis, A. & Kylafis, N.D. 1996, A&AS, 120, 171 (TMK96) 
Titarchuk, L., Mastichiadis, A. & Kylafis, N.D. 1997, ApJ, 487, 834 (TMK97) 
TuroUa, R., Zane, S., & Titarchuk, L. 2002, ApJ, 576, 349 

TuroUa, R., Zane, S., Zampieri, L. & NobiUi, L. 1996, MNRAS, 283, 881 (TZZN) 

Wood, K.S., Titarchuk, L., Ray, RS., Wolff, M.T., Lovellette, M.N., & Bandyopadhyay, 
R.M. 2001, ApJ, 563, 246 

Zel'dovich, Ya.B. &; Raizer, Yu. P. 2002, Physics of shock Waves and High-Temperature 
Hydrodynamics Phenomena (New York: Dover) 

A. Orthogonality of the eigenfinctions Rk{T) and analytical evaluation of the 



Wc remind the reader the procedure for proving the orthogonality of eigenfunctions Rk 
related to different eigenvalues A| of the boundary problem (27-28). One can check that 
equation (28) can be written in the self-adjoint form 



square norm Hf. 




(Al) 



where 



p(T) = [(/3-l)/C„,/3]rXT), 



(A2) 



and 



^(r) = r^i-^mP-^)ei. 



(A3) 



This preprint was prepared with the AAS macros v5.0. 
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Let us suppose that and are eigenvalues corresponding to the eigenfunctions Rn{T) 
and i?^ respectively. Then we use equation (Al) to write 



d 

d^ 

d 



p{t) 



Rn 



dRn 

~d^ 

dRm 

~d^ 



+ Xt^{T)Rn } = 0, 



+ Xltp{T)Rm > = 0. 



Subtracting the second from the first yields, upon integration, 

/•TO 

(A^-A^)/ i^{T)R^{T)R,r^{T)dT ^ p{t) 

Jo 



, , dR^ T-> / \ dR, 

Rn{T)—, Rm{r) 



dr 



dT 



To 



(A4) 
(A5) 

(A6) 



Because of homogeneous boundary conditions (28) the right-hand side of equation (A6) 
vanishes, and consequently we conclude that 



iyRni Rm) 



i>{T)Rr,{T)Rm{T)dT = 0, if A^ A^ 



(A7) 



This result estabhshes the orthogonahty of the family of eigenfunctions Rkir). Ordinarly, one 
must resort to numerical integration in order to compute the square norm — \ \Rk\\^ oi the 
eigenfunctions Rk{T). Suppose that A^ is the eigenvalue corresponding to the eigenfunction 
Rn{T), and A is an arbitrary value associated with the general solution R{X^,t) given by 
equation (27). Then, in analogy with equation (A6), we can show that 



Pir)[R'-t-Ruf] 

^{T)R{\T)Rk{T)dT ^ — 



To 







(A^ - AD 



(AS) 



The integral that we seek is obtained in limit A^ — > A^ However, in this limit, both the 
numerator and denominator in equation (A8) vanish, and therefore must use L'Hopital's rule 
to obtain 



Hn — \ \Rk 



lim p{t 



^ = lim 

A^-^I Jo 

OR dRk 



TO 



k 



ij{r)R{X,r)Rk{r)dr 

^ d^R ^ "° 
— Rk 



(A9) 



_d{X^) dr "dTd{X^) 

Because Rkij) goes to zero when r — > and r = Tq and dR^/dr — >• when r — >• (see Eq. 
28) our expression for reduces to 

dR{Xlr^)dRk{T^) 



^'==^^"°^- a(A^) dr ■ 
Thus there is no need to employ numerical integration to determine the values of 



(AlO) 
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Using the differential property of the confluent function $(a, b, x) (Abramowitz & Stegan 
1070, Eq. 13.4.8) and the condition (29) we find that 

= -3m„„,ro^/(^-^) + 2 + A^/?/3, /? + 3, -3m.„,(/3 - l)r,]. (All) 
The derivative 

a^^_7o) ^ ^^^+V(/?-i)(^/3)$. + 1 + ^2^/3^ ^ + 2, -SrhoutW - l)ro] (A12) 
can be calculated using the function '^{z) defined as '^{z) — r'{z)/r{z). In fact, 

where (a)„ = a(a + l)(a + 2). ..(a + n — 1), (a)o = 1 and 

(a); = (a)„[^(a + n)-^(a)]. (A14) 

One can calculate (a)'^ using the expansion of ^—function for large value arguments z ^ 1 
(Stegan & Abramowitz 1970, Eq. 6.3.18) 

*(z) ~ In z - l/2z - l/12z^ + l/120z^ - l/252z^ + ... (A15) 
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Effective optical depth T^ffg 



Fig. 1. — Plot of the photon index F = Fi vs. effective optical depth Teff^ — 3(/9 — l)7houtTo 
for various values of 
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1 1 1 r 




T 1 1 r 



T 1 1 r 



1 2 3 

DImensIonless time, u=ayNQct 



Fig. 2. — Plot of 'P{u) {'P{u) oc [r^/*^^~''^9P/9r[r, lijl^^o} for three monochromatic source 
distributions: /(r) = exp (rjC,), for r] = 1.5, 2.5 and distributed according to the first eigen- 
function -Ri(t) (see text). Plot for r/ = is not displayed here, because it is concentrated 
very close to u = in the presented scale. Calculations are made for (3 = 2, rhout = 0.8 and 
To = 2. 
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0.5 

Frequency vjv^ 

Fig. 3. — Plot of emergent photon redshifted line spectra. There are four spectra for 
monochromatic source distributions, /(r) = exp(77^), for = 0, 1.5, 2.5 and distributed 
according to the first eigenfunction R\{j) (see text). All these four spectra are related to 
through transformation (22) (see text). Calculations are made for /3 = 2, mo«t = 0.8 
and To = 2. The fifth photon spectrum (magenta) is for the asymptotic case of tq ^ 1 (see 
formula 50). The monochromatic 5— function source is located at ^ = ^* = 2, and /? = 2. 



